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o\ : 

In the context of inflationary scenarios, the observed large angle anisotropy of the Cosmic Mi- 
crowave Background (CMB) temperature is believed to probe the primordial metric perturbations 
from inflation. Although the perturbations from inflation are expected to be gaussian random 
fields, there remains the possibility that nonlinear processes at later epochs induce "secondary" 
■ non-gaussian features in the corresponding CMB anisotropy maps. The non-gaussianity induced by 

' nonlinear gravitational instability of scalar (density) perturbations has been investigated in existing 

literature. In this paper, we highlight another source of non-gaussianity arising out of higher order 
scattering of CMB photons off the metric perturbations. We provide a simple and elegant formal- 
ism for deriving the CMB temperature fluctuations arising due to the Sachs- Wolfe effect beyond the 
linear order. In particular, we derive the expression for the second order CMB temperature fluctu- 
ations. The multiple scattering effect pointed out in this paper leads to the possibility that tensor 
metric perturbation, i.e., gravity waves (GW) which do not exhibit gravitational instability can still 
contribute to the skewness in the CMB anisotropy maps. We find that in a flat Q, — 1 universe, the 
skewness in CMB contributed by gravity waves via multiple scattering effect is comparable to that 
, from the gravitational instability of scalar perturbations for equal contribution of the gravity waves 

and scalar perturbations to the total rms CMB anisotropy. The secondary skewness is found to be 
smaller than the cosmic variance leading to the conclusion that inflationary scenarios do predict that 
the observed CMB anisotropy should be statistically consistent with a gaussian random distribution. 
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I. INTRODUCTION 

X" 

Since its discovery by Penzias and Wilson Q, the Cosmic Microwave Background (CMB) has proved to be an 
extremely significant observational guide in our quest towards understanding the universe. The detection of tiny 
anisotropies in the CMB by the COBE - DMR group Q was an important milestone in the study of the universe and 
the understanding the large structure that we see around us. The COBE detection has opened up a fresh avenue of 
investigation and has been followed by a host of new developments both on the observational and theoretical fronts || . 

The idea of incorporating an inflationary phase in the early universe Q] has gained wide acceptance in the last decade 
and is perhaps the most prevalent scenario within which one attempts to understand the universe. It was realised soon 
after the notion of an inflationary scenario was put forward that besides resolving some long standing problems of the 
Big-Bang model of cosmology, inflationary models also predict the form of the power spectrum of the primordial scalar 
metric fluctuations which could seed the formation of the large scale structures observed in the present universe || . In 
fact, both gravity waves (GW), i.e., tensor metric fluctuations ||, as well as adiabatic density perturbations (related 
to the scalar metric fluctuations) arise as natural consequences of the inflationary scenario, due to the superadiabatic 
amplification of zero-point quantum fluctuations occurring during inflation. As gravity waves and scalar density 
perturbations enter the horizon, they induce distortions in the cosmic microwave background (CMB) through the 
Sachs- Wolfe effect 0. The relative contribution of the gravity waves and the adiabatic perturbations is linked to 
the specific model of inflation |J . The spectral index of the power spectrum of initial perturbations can be inferred 
directly from the CMB anisotropy measurements. 
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The measurements of the power spectrum of CMB temperature fluctuations have till date been found to be consistent 
with inflationary scenarios of the early universe. In particular, the spectral index inferred from the COBE - 4year data 
(see for eg. [|l0| ) is consistent with the near scale- invariant spectrum of fluctuations generically predicted by inflationary 
scenarios. Another generic prediction of inflation is that the metric perturbations generated are gaussian random 
fields. At the linear order, the CMB anisotropy produced would reflect the gaussian nature of initial perturbations. 
However, one cannot apriori deny the possibility that non-linear corrections to the growth of perturbations and to the 
Sachs- Wolfe effect could induce some non- gaussian features in the CMB even for gaussian initial metric perturbations. 

Non-zero skewness is a definite signature of non-gaussianity in a distribution. For gaussian initial perturbations 
the skewness in the CMB appears only beyond the linear approximation. At the leading order the skewness arises 
from two distinct effects. The first effect is that initially gaussian metric perturbations become non-gaussian when 
the lowest order nonlinearity becomes important in the course of their evolution due to gravitational instability. The 
non-linear component of the metric perturbations are non-gaussian and introduce non-gaussian anisotropies in the 
CMB through a linear Sachs- Wolfe relation at the corresponding order. The second effect is that the gaussian metric 
perturbations introduce non-gaussian anisotropies in the CMB due to second order (double) scattering of the photon 
off the linear order metric perturbations. This arises from the second order terms in the Sachs- Wolfe relation. All 
previous discussions of skewness in the CMB, have been limited to the the estimation of only the first effect, i.e., 
nonlinearity (and consequent non-gaussianity) due to gravitational instability The tensor component of 

metric perturbations (GW) does not exhibit gravitational instability, consequently the possibility of non-gaussianity 
in the CMB caused by the GW background has been entirely ignored in previous literature. 

In this paper, as an illustration of the second effect, we calculate the CMB skewness produced by a gaussian 
stochastic linear gravity wave background generated during inflation. In the context of f2 = 1, flat FRW models, the 
magnitude of the effect considered here appears to be comparable to the corresponding estimate of CMB skewness 
arising from the gravitational instability of scalar metric perturbations p2| . 

In §0, we outline the basic formalism involved in estimating the skewness in the CMB and present a very general 
approach for obtaining higher order corrections to the Sachs- Wolfe effect for a general cosmological perturbation. In 
the following section (£ III ) we estimate the skewness in the CMB anisotropy that would arise from a inflationary 
gravity wave background for a range of values of the spectral index. 



II. FORMALISM 



In this section, we outline the basic approach and present the derivations of results used in our calculation. The first 
part of the section contains a brief discussion of the perturbative approach used in estimating non-gaussianity in the 
CMB anisotropy. The second part gives a compact derivation of the CMB anisotropy arising through the Sachs- Wolfe 
effect upto second order in the primordial metric fluctuations. 



A. Non-gaussianity and Nonlinearity in the CMB anisotropy 

It is possible to address non-linear effects in the CMB within a perturbative framework by expanding the temperature 
fluctuations, AT/T, in orders corresponding to the powers of the initial metric perturbation as : 



AT _/AT\ (1) /AT\ (2) / AT\ 

— ~ y~ ) + v~/ + \~) 



(3) 

. (1) 



Given that the initial metric perturbations from inflation are linear and gaussian, any non-gaussian feature in the 
CMB maps can only arise from the higher order temperature anisotropy such as (AT/T)( 2 ). We shall call this higher 
order effect — secondary non-gaussianity. (The term "secondary" is used to denote the effects which take place after 
recombination. The effects prior to recombination are "primary".) 

A non-vanishing skewness is a definite signature of non-gaussianity in a distribution. At the linear order, the 
mean CMB skewness c| 3) (0) = ((ATW/T) 3 } = where () denotes averaging with respect to different realizations of 
stochastic space-time metric perturbations of the FRW cosmological model which produce AT/T. Substituting the 
expansion (|y) into the expression C3(0) = ((AT/T) 3 ), it is clear that the leading order (in powers of the initial linear 
metric perturbations) contribution is at the fourth order, 

(!) / AT\ (!) / AT\ ( 2 ) 
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(In this paper, we deal only with C3 (0) , and hence the superscript denoting the order has been dropped in the rest of 
the text) . It is clear from the above expression that the mean skewness depends not only on the mag nitude of AT^ /T 
but also on the extent of correlation of this term with the linear order terms, AT^/T. For example, in the case of 
scalar perturbations the second order terms which arise due to the non-linear evolution of density perturbations grows 
linearly with the expansion of the universe and can attain values AT^/T « 0.1 AT^/T at late times |J|. However, 
in a flat fi = 1 universe, the linear order term contributes only close to the surface of recombination (jj rj ICC ) 
and the second order term attains its largest value only at late times (77 sa 770). Consequently in the final result for 
the mean skewness, the decay of correlation between the linear and the second order term over this large physical 
separation (« rj ) along the line of sight attenuates the effect of the growth of the second order term, leading to a very 
modest value for C 3 (0) @. It was also pointed out in the same paper that the mean skewness which arises due to 
weakly non-linear density perturbations is expected to be somewhat larger in models where the linear gravitational 
potential changes at late times leading to a significant linear order integrated Sachs- Wolfe contribution at late times 
(eg. CDM+A, n + 1 models). 

Even in a flat = 1 universe, contribution to linear order integrated Sachs- Wolfe effect comes from inflationary 
tensor perturbations (gravity waves). Consequently, one expects that the correlation between the linear and second 
order terms is not attenuated in this case leading to larger values of 6*3(0). The second order AT*- 2 ' /T in the case 
of gravity waves comes only from double scattering since GW do not exhibit any gravitational instability. Scalar 
perturbations also give rise to second order anisotropy through double scattering. However, for flat, f2 = 1 models the 
contribution to the mean skewness is expected to be even smaller than that from gravitational instability considered 
in Jl^] . This can be seen from the fact that AT^ /T from double scattering too has contributions only at late times 
implying attenuated correlation with the linear term contribution close to the surface of recombination. 

B. Second order CMB anisotropy from the Sachs- Wolfe effect 

In a perfectly isotropic universe the CMB would have the same temperature in all directions on the celestial sphere. 
If, however, the cosmological metric is perturbed, the temperature observed today fluctuates over the celestial sphere. 

The dominant contribution at large angular scales (8 > 1°) to the observed temperature fluctuations comes from 
the change in the frequency of any CMB photon as it travels from the surface of last scattering to us []. In the case of 
an isotropic universe, the overall increase in the scale factor a(j]o) / a(rj rcc ) redshifts the entire Planckian distribution 
of photons leading to a Planckian distribution at a lower temperature given T roc /T = a(?7o)/a(?7 rec ). The presence of 
the perturbations h a b(r),x), produces an additional change in the frequency and direction (momentum) of a photon 
as it moves in and out of the fluctuating metric perturbations. 

We consider the trajectory of a photon (or ray) in a perturbed flat FRW universe and work in a synchronous 
coordinate system where the line element is of the form 

ds 2 = a(r]) 2 [~d V 2 + (S ab + h ab (77, x)) dx a dx b ] . (3) 

In the above, h a {> = h& ?(<~ e) + g2 ) ^ s the metric perturbation, and e < 1 is a small number characterizing the 

amplitude of deviations from the unperturbed background FRW universe. 

The photon trajectory can be obtained by perturbatively solving the eikonal equation for the phase S(x,rj). The 
eikonal equation for a photon propagating in a spacetime with metric, g^, is 

dSdS f» = 0, (4) 



dx^ dx u 



(analogous to the Hamilton- Jacobi equation for a massive particle). The frequency and the direction of the photon 
can be obtained from the phase S(x, rj) using 

1 d 

u{x,rj) = --^-—S{x,rf), k a {x,rf) = W a S[x,rf) . (5) 

Retaining terms to order e 2 when inverting the metric given by the line element in equation. (|J), the eikonal 
equation (||) becomes 



1 We assume instantaneous recombination which is an excellent approximation to standard recombination for calculating the 
CMB anisotropy at large angular scales (9 > 1°). 
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(6) 



where we use the background spatial metric 5 a b to raise and lower the spatial indices. The phase S(x, rf) can be 
expressed in a perturbative expansion in powers of e, as 

S = S^+S^(~e) + SW(~e 2 ) . (7) 

and substituting the above in equation (||) we obtain the zeroth order equation 

^-S^Y + V a S^V a S^ =0, (8) 
drj J 

the first order equation 

_ |. 5 (o)|_5(i) + Va ^°)v a 5« = ift«° 6 V a 5Wv 6 5W , (9) 
and the second order equation 

_ S (.0)B S W + V a S^V a S^ = -( |-S (1) V - -V S«V°S« + h^VaSWViSW 
or) or] 2 \or] J 2 



hm a b _ h(1) ac h(1)b 



V a S (0) V h S (0) . (10) 



The solutions to these equations correspond to a family of rays (or null geodesies) in a perturbed FRW universe. 
The solution to the zeroth order solution is 

S^(x,r,)^k a x a -r, + C . (11) 

This corresponds to a family of trajectories for the photons with frequency l/a[rj) moving in the k direction. In 
figure 0, we show some of these trajectories for the case where k is along the x axis (the y direction has not been 
shown). The observer O sits at the origin of the spatial coordinate system and measures the frequency of the photons 
at the present epoch 770. Using A as a parameter along the rays, this family of rays can be expressed as 

x(A) = kA + X (12) 

and 

r]{X)=X + r hn (13) 

with the condition that X is perpendicular to k. Here different values of (r) in , X) correspond to different rays and 
(?7o , 0) corresponds to the ray that the observer O sees at present traveling in the k direction . 

Substituting the first order solution in equation (0), we rewrite the LHS as a derivative along the zeroth order 
trajectory, i.e., with respect to the parameter A as 

J^«(z(A), 77(A)) = lhW ab (x(\), v (\))k a k b (14) 
which can be easily integrated to obtain the first order solution, 

S^i^xiXr)) = \ ^ h$(r,(\),x(\))k a k b d\ . (15) 



The frequency of the photon is given by the time derivative of given by equation (15). At the point A (refer 
to figure |l|) the time derivative is the A77 — > limit of 

(16) 

where A77 is the difference in r\ between the points B and A. 
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The point A lies on the trajectory x = kA, 77 = A + 770 and it corresponds to a value Ai for the parameter A. The 
point B lies on a different trajectory x = kA, 77 = A + r]o + A77 and it corresponds to the same value of the parameter 
i.e., Ai. We then have 
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S(B)-S(A) = - I ^ 1 fc ) ( ?? (A) + Ar ? , a; (A))-C(r ? (A), a; (A)) 



k a k b d\ 



(17) 



which leads to 

l-SWMAO^Ai)) = \£ ^ b \x)k a k b dX . 

Similarly, in the direction parallel to k (the x-axis in fig.|l|), the spatial derivative is given by limit 

S(C) - S(A) 
Ax 



(18) 



(19) 



where Ax is the difference in x between the points C and A. C lies on a trajectory x = kA, r\ = A + ?/o — Ax and it 
corresponds to the value of the parameter Ai + Ax. We obtain 



S(C) - S(A) = 



1 



(77(A) - Ax,x(A))fc Q fc"dA - - / h% (77(A), x(X))k a k b dX 



(20) 



and the spatial derivative in the direction parallel to k reads 



V c S^(r,(X 1 ),x(X 1 ))='f 



h ( ^ b \x(X 1 ), v (X 1 ))k a k b 



Xl ^(\)k«k b d\ 



(21) 



In the direction perpendicular to k we have 



V c 5( 1 )(x(A 1 ),t ? (Ai)) = \ V c h^(X)k a k b dX 
* Jo 



(22) 



Substituting the derivatives of in the second order equation (10) and using Vt a = V a — k a k b V b to denote the 
spatial derivative in the direction perpendicular to k we obtain 



■s" -'(.-•. //! = / dX i ^ ab {X)k a k b hW cd {\)k c k d - -^ ab {X)k a k b 



dud 



drj 



h {1 \ lj {X)k i k : >dX 



hW{\) - h^ u {x)h^\{x)\ + hi^ a \\)k a I v_ Lb h^ lJ {\)k i wd\ 



I (\^ l \ ] (X)¥yd\ \ [ \ [ V ± a h^ lm (X')k l 



2 Jo 



2 Jo 



l k m dX 



(23) 



The perturbative solution for <S*(x, 77) can now be used to obtain the frequency of a photon at any point along 
its trajectory. Although we are interested in a photon that leaves the last scattering surface at 77^ and reaches the 
observer at 770 , it is convenient to consider a photon traveling backwards in time from the observer to the last scattering 
surface. Keeping terms upto order e 2 we find that the photon that is observed at the frequency u> = l/a(i] ) left the 
last scattering surface with the frequency, 



a(77 e 



(24) 



Invert the above relation we find that a photon which left the last scattering surface with the frequency lo = 1/ a(j] e ) 
will have frequency 



1 



1 + ^(^e) + l^^e) + {^{x^f 



(25) 
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when it reaches the observer. This relates the observed frequency to the emitted frequency and the metric perturba- 
tions. Using equations (|24| ) and (^5| ) we obtain expressions for the fractional change in the frequency of the observed 
photon relative to the frequency that would be observed if the universe were unperturbed. Since the CMB photons 
have a Planckian distribution, (frequency independent) fractional changes in frequency translates to fluctuations in 
the temperature characterizing the distribution. At the linear order we recover the familiar (linear order) Sachs Wolfe 
effect 



AT' 1 
T 



(26) 



At the second order the expression for fractional change in CMB temperature reads 



AT(2 
T 



^ dX \ jh^ ab (X)k a k b -^h^ cd (X)k c k d 

A 



ld_ 
Adr)' 



Va dT ) 



h {1] ab {X)k a k b [ ^-h^ l0 {X)k l ¥dX --h {1 \ b {\)k a k b [ ^h^ijiX^VWdX' 



drj 2 




t lL h W a (\)k a J V Lb hW lJ (X )k t k^dX +h^ a (X)k a j ^-V La hW lJ (X )k'k^dX 
^V La h^ tJ (x')k t k^dx) Uj X V L b h^ lm ,(X )k l k m dX 



i rVo f) 



(27) 



The expression for AT^ /T has been independently obtained by other methods p7| , ^6| . 

The AT/T calculated here is the fractional increment (or decrement) with respect to the CMB temperature that 
would be observed in a homogeneous and isotropic universe. In practice an observer would measure the fluctuations in 
the CMB temperature from different directions in the sky and in order to relate our calculation with what is observed 
we have to remove any monopole component in the angular distribution of the calculated AT/T. In addition, the 
dipole component of AT/T is usually interpreted as arising from the observers peculiar motion and this is removed 
before relating the observed AT/T to the primordial fluctuations. 

In most previous work on second order CMB anisotropy, only the term depending on arising from nonlinear 
gravitational instability of scalar density perturbations in the above expression for AT' 2 ' jT has been considered. The 
nongaussianity from the term in AT^> /T for fl = 1 models has been considered in reference [|l2|,|14||. The 



term in AT^/T for < 1 models has been calculated in reference 18 to calculate the fourth order corrections to 
the power spectrum of CMB anisotropies. 



III. SKEWNESS FROM GRAVITY WAVE 



Gravitational waves can be decomposed into two polarization states denoted by + and x , and the metric pertur- 
bation corresponding to a gravitational wave with wavelength ^ traveling in the z direction can be written as 

hi b ( V , X )=e^h q (r ] )[a+e+ + aje*] ( 28 ) 

q2 

where the quantity A^(q)/q^ is the power spectrum of the GW perturbations. The temporal evolution of the modes 
of the gravity waves in a fi = 1, matter dominated FRW universe is given by 

h q (v) = —Ji(qv) , (29) 
qi] 



G 



where j\ denotes the spherical Bessel function of order one. In eqn. (|2q), the a^s are the complex amplitudes of the 
two polarization states (henceforth denoted by the superscript a) and e" b are the usual basis tensors for the traceless 
metric perturbations. 

A gravity wave in any arbitrary direction n can be obtained by rotating the coordinate system (or equivalently the 
wave traveling in the z direction) so that the wave travels in the required direction. Denoting the rotation by R(n) a b, 
we can express a gravitational wave traveling in the ii direction as h™ b = R c a (n)R b {n)h z ab . 

Thus we can write metric perturbations corresponding to an isotropic stochastic background GW as 



h ab (77, x) = 



: (n)R db {n)[a+e+ d + a*e c * ]A G {q) 



(30) 



In the above, a" are independent random variables which satisfy the relation 



<<aV >=<5 J (q-q)<S' 



a/3 



(31) 



where the angular bracket denotes an ensemble average over different realizations of the stochastic GW background. 

We substitute the expression (|3(]) for the GW metric perturbation in equations (|26| ) and ( p7|) to calculate the 
temperature fluctuation in the CMB at the first and second orders, respectively. At the linear order, the AT/T 
produced by gravity waves has no monopole or dipole component, but at the second order there is a monopole 
component which has to be removed. To keep the algebra at a manageable level we use the ensemble average of 
AT^/T instead of the angular average of AT^/T to estimate the monopole component and we subtract this 
contribution from the second order temperature fluctuation predicted by equation ( |2"7| ) . We next use this to calculate 
the skewness. 

At the leading order the skewness is given by 



C 3 (0) = 



/AT\ 3 \ _ _ / /AT (1) 



V T 



/AT(2) 



(32) 



and substituting for expressions for ATW/T and AT^ /T we obtain 

3 5 f0~c ^ fVrec _ fVtec ^ fOO ^ ,O0 dq 



C 3 (0) 



2 5 7T 10 



j2(g2A 3 ) j2(g2A 2 ) 



A 3 
3 5 

2 5 7T 10 



A 2 



dX 2 

3?r 
T 



dX 3 



no 



qi Jo 92 



Aa{qi)Ah(q2) 



ii(9iA 3 ) J2(9iAi 



9i A 3 



Ai 



J 3 (9i(A 3 - Ai)) J 3 (gi(A 3 - A 2 )) _ j2(gi(A 3 - Ai)) j 2 (gi(A 3 - A 2 )) 
(9i(A 3 -Ai)) 2 (gi(A 3 -A 2 )) 2 (^(As-Ax)) 2 (gi(A 3 -A 2 )) 2 

j2(qi(h - Ai)) j 2 (92(A 4 



' r5 ° f°° d 92 .2 , n A 2 I 

— \ — ^g(9i)^g(92 



dA 3 / dXi± 
I a 3 ' > 



9i 



A2)) 



j2(qiXi) j 2 (q2X 2 ) 



Ai 



A 2 



32(qih) h{q2^A) ji(giA 3 ) g 2 . 4 



A 4 



<7iA 3 A 4 



92 

j 2 (g2A 4 ) 



(9i(A 3 -A!)) 2 fe(A 4 -A 2 )) 5 



A 4 



ji (92 A 4 



+ 



2 5 7T 10 



dXi 



dX 2 



dq ,2 ,^j2(q(\ 



... ~ Ai)) j 2 (gAi) j 2 (giA 2 ) 

q Ag{Q > (q(X 2 Ai)) 2 



Ai 



A 2 



6C 2 (0), 



(33) 



where j n and J n denote spherical and ordinary Bessel functions of order n, respectively. In the case of matter 
dominated f2o = 1 models with three relativistic species of neutrinos, 77 rcc = ?/o/49.6. The final results are however 
completley insensitive to the exact value of ry rcc . It is interesting to note that none of the terms in equation (|27]) which 
have spatial derivatives in direction perpendicular to k contribute to the skewness. 

We numerically evaluated the above expression to compute the skewness for different spectral indices of GW power 
spectrum. A dimensionless skewness parameter S3 can be constructed by dividing the skewness, C 3 (0) by the square 
of the variance, C 2 (0). The variance of AT/T that arises due to relic gravity waves from inflation is given by 



C 2 (0) = 



2V 



dX x 



dX 2 



dq . 2 M hWte 
9 (9(A 2 - 



A i)) j2(gAi) j 2 (giA 2 ) 



Ai)) 2 Ai 



A, 



(34) 



The power spectrum of the initial gravitational wave perturbation is assumed to be a scale-free power law, A^q{U) — 
A k nT , where tit = corresponds to a scale-invariant spectrum. Power law models of inflation produce gravity waves 
with riT < 0. We compute the CMB skewness for a broad range in tit (—0.5 < n-r < 0). We find that the value of 
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S3 varies from —1.6 to —3.1 for n<r varying from the scale invariant spectrum to Ut = —0.5. It is interesting to note 
that the skewness arising from gravitational instability of scalar perturbations is comparable for the same range of 
tilt ( S3 w —2.2 for scale invariant spectrum p2[ ). 

The corresponding observable quantity for the skewness, 63(0) is the all sky-average C^O) = 
(47T)- 1 / {AT{6,<p)/TY dVL of one particular realization of the random fluctuations. The value obtained by tak- 
ing an angular average over one sky would generally differ from the ensemble average over all realizations by the 
(cosmic) variance of the observed skewness, 6*3(0). The skewness originating due to any effect would have an ob- 
servable significance if the predicted signal stands above the cosmic variance. This is a fundamental limitation and 
a minimal requirement. In practice, a detectable signal has to stand above additional variances such as instrumental 
noise, finite beam width of antennas, incomplete sky coverage etc. 

The Cosmic variance can be expressed in terms of an angular integral over the cube of the two-point correlation 
function, Ci (9) pl| . Assuming a gaussian approximation for the two-point correlation function, we express 

C 2 {6) = C 2 (0) exp[- ^ c + 1} 02 ], C 2 (0) « 3 x lO" 5 (35) 

where the cut-off, l c , in the AT/T angular spectrum at large values of the spherical harmonic eigenvalue, I (l c w 
Vo/Vrec ~ 49 for GW). Using equation (p5|), the Cosmic variance, SS3, for the case of a CMB anisotropy arising from 
gravity waves is given by SS3 w l/(C2(0jl c ) « 670. The Cosmic variance for the scalar case is around 5 times smaller 
|L2[ . The Cosmic variance is larger as the power spectrum tilted away (reddened) from scale invariance for both scalar 
perturbations and gravity waves. It is clear that in principle the secondary skewness in the CMB for a CDM model 
(f2 = l,fli = 0.05 and Hq — bOkms / s / Mpc.) is unobservable since it is below the cosmic variance. 



IV. CONCLUSIONS 



In this work we study a possible source of secondary non-gaussianity in the CMB and reexamine whether initial 
gaussian metric perturbations (as expected from inflation) should lead to a gaussian CMB anisotropy. We point 
out that besides nonlinear gravitational instability, secondary non-gaussianity can be induced in the CMB maps due 
to multiple scattering of CMB photon off metric perturbations. We develop a general method for calculating CMB 
temperature fluctuations beyond the linear order and calculate the expression for CMB temperature fluctuations at 
the second order. As an example of skewness from multiple scattering, we estimate the skewness in the CMB that 
is expected to arise due to a relic inflationary gravity wave background. We show that the magnitude of the effect 
studied here is comparable to the ones considered in existing literature. However, the signal is much smaller than its 
cosmic variance and we conclude that gaussian initial perturbation from inflation leads to a gaussian statistics for the 
observed temperature fluctuations. 
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FIG. 1. A sketch of the unperturbed rays of CMB photons propagati ng: al ong the x axis is shown (the third spatial dimension 
is suppressed). The figure serves as a guide to the arguments used in STIB to compute the temporal and spatial derivatives of 
the phase S(rj,x). The figure clearly illustrates the distinction between spatial derivative taken along the direction of photon 
propagation and spatial derivative taken perpendicular to the direction of photon propagation. 
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